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ROTATIONAL  ELASTIC  DYNAMICS 


by 


J.  Baillieul  and  M.  Levi 

Aerospace  and  Department  of  Mathematics 

Mechanical  Engineering  Boston  University 

Boston  University  Boston,  MA  02215 

Boston,  MA  02215 


Abstract 

The  combined  dynamical  effects  of  elasticity  and  a  rotating  reference  frame  are  ex¬ 
plored  for  structures  in  a  zero  gravity  environment.  A  simple  yet  general  approach  to 
modeling  is  presented,  and  this  approach  is  applied  to  analyze  in  detail  the  dynamics  of 
a  specific  prototypical  structure.  Energy  dissipation  is  included  and  its  effects  are  studied 
in  detail  in  a  model  problem.  Bifurcations  and  stability  axe  analyzed  as  well. 
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1.  Introduction 


There  is  now  a  fairly  general  awareness  among  aerospace  engineers  that  the  dynamics 
of  the  complex  spacecraft  currently  in  production  and  on  the  drawing  boards  will  be 
greatly  influenced  by  continuum  mechanical  effects  such  as  elasticity.  Indeed  as  the  designs 
being  contemplated  increase  in  size  and  complexity  ([13])  the  dynamic  effects  of  flexible 
members  become  more  important,  and  recognizing  this,  many  researchers  over  the  last 
decade  have  focused  their  efforts  on  obtaining  new  methods  for  the  design,  analysis,  and 
control  of  flexible  mechanisms.  Space  does  not  permit  (nor  would  it  be  in  keeping  with 
the  main  purpose  of  this  report)  to  survey  the  vast  literature  on  flexible  space  structures: 
the  interested  reader  can  get  some  idea  of  research  activity  in  this  area  by  referring  to  any 
one  of  a  number  of  collections  of  papers  and  conference  proceedings,  6uch  as  [13]  or  the 
more  recent  volume  [11]. 

The  purpose  of  this  report  is  to  describe  recent  research  which  has  been  aimed  at 
developing  a  mathematical  theory  of  the  rotational  dynamics  of  complex  mechanical  sys¬ 
tems  which  include  articulated  and  elastic  components.  Our  objective  in  this  research  has 
been  to  carry  out  a  study  of  the  global  qualitative  dynamics  of  such  systems  in  sufficient 
depth  as  to  allow  predictions  regarding  the  stability  and  asymptotic  behavior  of  spacecraft 
due  to  a  variety  of  energy  dissipation  mechanisms  such  as  viscoelastic  material  damping 
of  vibrations  of  clastic  parts.  We  believe  that  historical  evidence  points  to  the  value  of 
developing  a  fairly  complete  global  asymptotic  stability  theory  of  this  type  since  there 
are  numerous  examples  of  missions  in  space  which  did  not  achieve  their  stated  objectives 
because  certain  long  term  mechanical  effects  were  never  adequately  taken  into  account  in 
the  mission  planning.  Explorer  I,  the  first  successfully  launched  American  satellite,  pro¬ 
vides  the  best  known  example  of  such  untoward  behavior.  Upon  achieving  earth  orbit, 
the  pencil-shaped  satellite  was  supposed  to  rotate  about  its  major  axis  of  inertia.  Before 
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it  had  completed  one  orbital  revolution,  however,  radio  signals  indicated  that  a  tumbling 
motion  had  developed  and  was  increasing  in  amplitude. 

While  explanations  of  Explorer’s  errant  behavior  have  been  offered  by  a  number  of 
researchers  (see  e.g.  references  to  this  problem  in  [7]),  we  axe  aware  of  no  attemps  to  obtain 
a  rigorous  mathematical  analysis  of  such  occurences,  with  the  exception  of  [3],  where  some 
of  the  results  of  the  present  report  were  announced.  We  also  mention  a  recent  paper  [8] 
which  deals  with  a  somewhat  different,  although  related  aspect  of  a  similar  problem  and 
several  more  recent  efforts  reported  in  [111  which  develop  results  similar  to  those  described 
in  Section  5  below. 

Our  report  is  organized  as  follows.  In  Section  2  we  derive  equations  describing  the 
rotational  dynamics  of  complex  structures.  These  include  the  general  effects  of  inertial 
forces  created  by  rotation  of  the  reference  frame.  Section  3  focuses  on  a  general  theoretical 
framework  for  Lagrangian  mechanics  with  damping.  In  Section  4,  a  simple  structure 
consisting  of  a  rigid  body  with  an  elastic  beam  appendage  is  studied,  and  we  present 
what  we  believe  is  the  simplest  reasonable  continuum  mechanical  model  of  such  a  system 
undergoing  three  degree  of  freedom  rotations.  The  asymptotic  steady  state  dynamics  for 
this  system  are  studied  in  Section  5,  and  in  Section  6  we  present  a  detailed  analysis  of 
the  qualitative  dynamics  of  a  closely  related  model  having  only  one  rotational  degree  of 
freedom.  A  second  part  of  the  report  treats  “Equilibrium  Mechanics  of  Rotating  Systems,” 
dealing  principally  with  rotating  kinematic  chains. 
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2.  The  rotational  dynamics  of  complex  structures. 


In  this  section  we  will  derive  equations  of  motion  for  a  class  of  structures  consisting 
of  elastic,  fluid  or  rigid  components.  While  these  equations  are  completely  general,  they 
are  most  useful  in  describing  any  structure  whose  configuration  is  conveniently  specified 
by  the  position  of  the  structure  with  respect  to  a  moving  coordinate  frame  together  with 
the  position  of  that  frame  in  space,  i.e.  relative  to  some  fixed  inertial  system. 


Figure  2.1:  The  body  frame  is  translated  and 
rotated  with  respect  to  the  inertial  frame. 

We  fix  an  orthonormal  basis  forming  an  inertial  frame  (“space  frame”),  and  choose  a 

“body”  coordinate  system  designated  by  a  set  of  orthonormal  vectors  gi,g2,g3*  Choice 
of  the  gi’s  will  depend  on  the  particular  problem  at  hand,  and  it  affects  the  simplicity 
of  resulting  equations.  In  section  4,  where  we  treat  an  example  of  a  rigid  body  with  an 
elastic  appendage,  we  affix  the  body  frame  g<  to  the  rigid  body,  although  other  choices  are 
possible  [6]. 

The  position  and  orientation  of  the  body  frame  (with  respect  to  the  chosen  inertial 
frame)  may  be  described  at  each  time  t  by  an  element  of  the  special  Euclidean  group, 
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SE( 3,R),  represented  by  a  4x4  matrix 

m-(r®  f) 

(see  figure  2.1),  where  Y (t)  G  S0(3)  is  an  orthogonal  matrix  describing  the  orientation  of 
the  body  frame  and  y  G  R3  is  the  position  of  the  origin  of  that  frame  in  space,  (cf.  [5]). 

Therefore,  if  a  point  P  is  given  by  the  vector  u  in  the  body  frame  and  by  vector  U  in 
the  space  frame,  then  u  and  U  are  related  via 


U  =  Xu, 


where  U  =  ( ^ )  and  u  =  ( “ ). 

The  positions  of  various  elements  of  our  flexible  structure  relative  to  the  body  frame 
will  be  described  by  a  vector  function  u(z,t)  =  u(zi ,  zj ,  zj ,  t)  denoting  the  position  at  time 
<  of  each  particle  whose  “unperturbed”  position  is  at  z.  Here  “unperturbed”  can  mean 
either  initial  or  undeformed,  depending  on  one’s  choice.  In  section  4,  z  will  denote  the 
neutral  position  of  a  particle  of  the  flexible  structure. 

In  summitry,  we  consider  systems  whose  configuration  space  is  given  by 
{g}  =  {(y,u,y)}  =  SE(Z)  xC  =  50(3)  x  R3  x  C,  where  C  =  {u(-,<)}  is  a  suitably  defined 
function  space  whose  elements  are  functions  u(z)  describing  the  configuration  of  the  body 
relative  to  the  body  frame. 

We  describe  now  the  kinematics  of  the  system  and  give  an  expression  for  the  kinetic 
energy.  The  evolution  of  the  matrix  X(t)  G  SE(Z)  can  be  described  by  a  differential 


where  ft(<)  is  the  skew-symmetric  matrix 


-u>3 

0 

u>i 
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of  angular  velocities  about  the  corresponding  body  axes  gj,  t  =  1,2,3. 

The  inertial  coordinates  U{z,i)  of  a  point  P  are  related  to  its  body  coordinates  u(z,t) 
via 

U(z,t)  =  Xu{z,t) 

and  the  corresponding  velocities  are  given  by 

jfU  =  Xu  +  Xut  = 

_  +  ut )  + 

where  Ut  denotes  the  partial  derivative  with  respect  to  t. 

The  kinetic  energy  of  the  body  is  then  given  by 

n«, «)  =  \  Jg  \\vrim  =  i  jf  ||y(n«  +  «,)  +  y\\2im,  (2.1) 

where  B  denotes  the  point  set  comprising  the  body  (at  time  t)  described  in  the  body 
coordinate  system,  and  dm  is  the  mass  distribution  in  this  coordinate  system. 

The  kinetic  energy  of  almost  any  rotating  structure  will  be  of  this  form,  as  it  does  not 
depend  on  the  constitutive  relations  governing  the  structure  itself.  More  refined  dynam¬ 
ical  models,  as  treated  in  subsequent  sections,  will  embody  structural  information  in  the 
expression  for  potential  energy 

V(,)  =  V(l>,y) 

(which  we  assume  to  be  independent  of  q).  Without  specifying  the  form  of  V  at  this  point, 
we  prove  the  main  theorem  of  this  section. 

Theorem  2.1.  Let  A  be  defined  by 

=  /  ((fiu  +  Ui)uT)a<fm, 

Jb 
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where  Ma  =  |(M  —  MT)  denotes  the  antisymmetric  part  of  a  matrix  M.  Equations  of 
motion  of  any  system  whose  configuration  space  is  {9}  =  {(Y,u,y)}  =  50(3)  x  C  x  R3  as 
above  are  given  by 


A(t)  +  1(1(1),  A(t)\  -  J  (YTjuT).im  =  T-  (Y-'Vy). 

utt  +  n\  +  flu  +  2flu(  +  YTy  =  F  —  ~~ 

ou 

f  8V 

I  F(fl2u  +  flu  +  2fiut  4-  u«)  +  ydm  —  f  — — 

Jb  oy 

Y  =  YS} 


(2.2) 

(2.3) 

(2.4) 

(2.4) ' 


where  Vy  is  the  matrix  of  partial  derivatives  and  is  the  Frechet  derivative  of  V 

with  respect  to  u. 


In  the  skew- symmetric  matrix 


(  ° 

-r3 

T2  \ 

T  = 

TS 

0 

-n 

\-r2 

Ti 

o,  > 

Ti  is  the  net  nonconservative  torque  applied  to  the  system  about  the  body  axis  g F  =  F(z ) 
is  the  distributed  nonconservative  force  density  acting  on  the  particle  positioned  at  u(z,t) 
expressed  in  the  body  coordinate  system,  and  f  is  the  net  nonconservative  exogenous  force. 

Remark:  In  subsequent  sections  nonconservative  forces  will  arise  due  to  viscoelastic 
damping. 

Before  giving  the  proof,  we  rewrite  these  equations  so  as  to  provide  a  clearer  picture 
of  the  physical  situation.  Let  u;  =  (uji,u;2,u>3)t  be  the  angular  velocity  of  the  body  frame 
expressed  in  that  frame ;  u>  is  related  to  the  angular  velocity  matrix  SI  =  Y~lY  as  follows: 


n  = 


o 

U>3 

— OJj 


-U>3 

0 

u>i 


U>2 

-U,J 

0 


,  w  = 


(U>J 

U>2 

U)3 
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Denote  by  S  the  operator  taking  a  skew-symmetric  matrix  fl  into  vector  w  :  SQ  =  w.  One 
easily  verifies  the  following: 

Lemma  2.1.  Given  any  pair  A,B  of  skew- symmetric  3x3  matrices  and  any  pair 
u,v  of  S-veciors,  the  following  identities  hold 

(i)  S([A,Bl)  =  S(.A)xS(B) 

(ii)  S(uv  r  —  vuT)  =  v  x  u 

(iii)  Au  =  5(A)  x  u 

Applying  the  operator  5  to  equation  (2.2)  and  using  Lemma  2.1  on  equations  (2.2)-(2.4) 
we  obtain 


Corollary  2.1  Equations  of  motion  (t.t)-(t.4)  are  equivalent  to 


a(t)  +  u>(<)  x  a(t)  4  /  u(z,t)  x  Y^y(t)dm  =  S(T  -  (Wy).)  (2.5) 

Jb 


,  .  ,  -  fV 

utt  -(-  w  x  (w  x  u)  +  w  x  u  -(-  2w  x  V  y  —  F  — — 

ou 


(2.6) 

(2.7) 


r  dv 

J  \Y(uj  x(wxa)+wxu  +  2wxui  +  Utt)  +  y]dm  —  f  -  — 

where  a  is  defined  by 

a(t)  =  /  ux(ui+wx  u)dm. 

Jb 

Written  in  this  way  equations  (2.5)-(2.7)  give  an  explicit  description  of  the  inertial 
forces  on  the  mechanical  system  viewed  in  the  moving  body  frame.  Introducing  the  deriva¬ 
tion  D  =  ji(-)  +  v  x  (•))  we  obtain  yet  another  rendering  of  these  equations. 

Corollary  2.2  Equation  of  motion  (t.t)-(t.4)  are  equivalent  to 


f  u  x  {D3u  +  Y~iy)dm  =  S{T  -  (Y^Vy),.) 
Jb 


(2.5)' 
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^  I  3 

Lq  |<^o 

1 

II 

(2.6)' 

i 

II 

J 

(2.7)' 

The  motions  of  any  complex  structure  undergoing  free  or  forced  rotation  are  described 
by  equations  (2.5)-(2.7).  This  formulation  is  thus  fairly  general,  and  it  cam  incorporate 
external  forces  and  torques  (due,  for  example,  to  gravitational  and  magnetic  fields)  and 
interna’  forces  (due,  say,  to  actuation  of  joints  or  the  constitutive  properties  of  the  mate¬ 
rial).  In  the  next  section  3  we  will  incorporate  dissipative  effects  in  this  formulation,  and  in 
section  4,  we  shall  develop  a  complete  dynamical  model  of  a  system  where  the  constitutive 
relations  are  those  of  a  simple  damped  beam. 


Proof  of  theorem  2.1:  The  proof  will  be  given  in  two  parts:  First,  the  treatment 
of  the  inclusion  Y  6  SO(3)  as  the  holonomic  constraint  onto  50(3)  as  a  submanifold 
of  GL{ 3),  with  an  appropriate  modification  of  the  Lagrange  equations  (Lemma  2.2),  and 
second,  the  utilization  of  the  left-invariance  of  (the  rigid  path  of)  kinetic,  energy  to  simplify 
the  resulting  equations  (Lemma  2.4). 

We  would  like  to  write  the  equations  of  motion  of  the  structure  in  the  Lagrangian 


form 

d  dT  dT  dV 

He}-  a7  =  ~a??  =  (: 

An  appropriate  modification  of  the  Lagrange  equations  is  needed,  however,  to  account  for 
the  fact  that  Y  is  constrained  to  the  submanifold  S0(3)  C  Gl( 3).  Such  a  modification 
is  unnecessary  if  the  Lagrangian  equations  are  expressed  in  terms  of  a  local  coordinate 
system  on  the  constraint  manifold-cf.  [1],  page  77.  This  approach  ignores  the  symmetry 
in  our  system  of  equations,  however.  The  following  Lemmas  2.2  and  2.4  are  key  to  our 
proof  of  Theorem  2.1 
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Lemma  2.2  Any  extremal  q(t)  of  the  action  J  L(q,q)dt  with  q  constrained  to  a  sub¬ 
manifold  Mq  of  a  Riemannian  manifold  M  satisfies  the  differential  equation 


(±9L_  9L\ 

V dt  dq  dq  )  ~  °’ 


where  ng  is  the  normal  projection  from  TqM  onto  TqMo. 

We  omit  a  straightforward  proof.  Corresponding  to  the  three  components  of  q 
(Y,u,y)  we  obtain  three  components  for  the  equations  of  motion: 


'  d  dr 

9T 

av\  n 

(2.8) 

dt  dY 

~dY  + 

dYj 

d  ST 

ST  _ 

SV 

(2.9) 

dt  Su 

Su 

Su' 

d  8T 

8T  _ 

dV 

(2.10) 

dt  dy 

dy 

dy' 

where  Py  is  the  orthogonal  projection  from  Ty  GZ(3)  onto  TySO( 3)  in  the  trace  norm 
( A,B )  =  tr  AtB.  Here  y  are  the  derivatives  with  respect  to  the  standard  Rie¬ 

mannian  structure  on  TGL{ 3);  they  can  be  represented  as  matrices  of  partial  derivatives: 
a?  =  ^  =  ( efiF7)’  fu  denotes  the  Frechet  derivative  of  T  with  respect  to  the 

distributed  parameter  u.  The  expression  for  the  projection  Py  is  provided  by 


Lemma  2.3.  For  any  Y  £  50(3)  and  A  C  TyCl( 3) 

PyA  =  Y{Y~'A)ay 
where  Xa  =  \{X  -  XT) 

mt 

is  the  anti- symmetric  part  of  X. 


(2.11) 


Proof.  This  follows  from  decomposing  the  Lie  algebra  of  <j/(3)  as  the  orthogonal 
direct  sum  of  symmetric  and  skew  symmetric  matrices.  (gl( 3)  =  so(3)  ©  | 


Lemma  2.4  Suppose  T  is  a  left-invariant  function  on  TGl(n),  i.e.  suppose  there  is 
a  function  K  defined  on  gl(n )  (=  the  Lie  algebra  of  Gl(n)  =  space  of  real  n  x  n  matrices ) 


such  that  T(Y,Y)  =  K((l)  where  Q  =  Y^Y.  Then  for  (Y,Y)  e  TSO(n ) 

where  M  =  (§^)o>  and  ^  is  the  derivative  of  K  with  respect  to  Q  evaluated  at  (l  =  Y~1Y. 

Remark.  There  is  an  orthogonal  direct  sum  decomposition:  5/(71)  =  so(ti)  0  «o(n)x , 
where  so(n )  =  the  Lie  algebra  of  n  x  n  skew-symmetric  met  rices,  and  where  orthogonality 
is  defined  in  terms  of  the  trace  form  inner  product  [A,B\  =  tr  ABT  defined  on  5/(71).  If  the 
function  K  appearing  in  the  statement  of  Lemma  2.4  can  be  decomposed  as  K  =  Ki  + 
where  K{  depends  only  on  the  i-th  component  (t  =  1,2)  in  this  orthogonal  direct  sum, 
then  M  =  hTj(17),  where  by  K[  we  mean  the  derivative  of  K\  with  respect  to  the  natural 
differentiable  structure  on  so(n)  defined  in  terms  of  the  Killing  form. 


Proof.  A  JT  and  §y  may  be  thought  of  as  elements  in  the  cotangent  bundle  T*Gl(n). 
Making  the  usual  identifications,  the  standard  Riemannian  structure  on  Gl(n)  may  be 
prescribed  explictly  in  terms  of  the  trace  form,  and  we  may  write  §y  and  ^  all  as 


7i  x  n  matrices 


£0T 

dtdY 


'  d  dK 
dt^dSl 


)-nT 


dK' 

an 


dT 

dY 


~dK 


For  (Y,Y)  €  T  SO(n),  we  have  (y-1)T  =  Y  and  flT  =  -0,  and  the  result  follows  from 


Lemma  2.3. 


I 


The  proof  of  Theorem  2.1  proceeds  as  follows.  Using  left-invariance  of  the  first  term  in 
the  expression  for  kinetic  energy  (cf.  (2.1)) 

T=  \  [  ||ftu  +  u,||2  +  2(ftu  -)-  uuY-'y)  +  ||y||2  dm, 

2  J  B 
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we  obtain  equation  (2.2)  as  a  consequence  of  Lemma  2.4.  The  remaining  equations  (2.3) 
and  (2.4)  follow  by  direct  computation  from  equation  (2.9)  and  (2.10).  We  omit  the  details. 
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3.  Lagrangian  Mechanics  with  Damping 

A  major  advantage  of  Lagrangian  versus  Newtonian  mechanics  is  the  invariance  of 
Lagrange’s  equations  with  respect  to  coordinate  changes;  it  is  this  invariance  that  facilitates 
significantly  the  derivation  of  our  equations  of  a  motion.  It  is  thus  desirable  to  have  the 
extension  to  the  dissipative  case.  Such  a  modification  is  described  in  [9],  [12];  we  reproduce 
it  here  in  a  slightly  more  general  form. 

Let  D(q,  q)  be  the  so-called  dissipation  function  defined  as  follows: 

qDg  =  rate  of  dissipation  of  energy  per  second. 


(One  can  think  of  Dq  as  the  generalized  dissipation  force,  and  q  =  velocity.  The  above 
simply  says:  velocity  •  force  —  power). 

Let  L(q,  q)  be  the  Lagrangian  of  the  system.  The  equations  governing  the  system  are 

jtL4  -  Lq  +  Di  =  0.  (3.1) 

Remark  3.1.  If  D  is  quadratic  in  q  (as  will  be  the  case  in  the  application  presented 
in  the  next  section),  then 

D  =  \qDq  =  ^(rate  of  dissipation). 

Equations  (3.1)  are  consistent  with  the  definition  of  D ,  as  shown  by 
Theorem  3.1.  If  E(q,q)  is  the  total  energy  of  the  system,  then 


Proof.  E  is  given  by  the  Legendre  transform*  in  q  of  L: 

rr  dL  r 

E  -  q—  -  L. 

8q 


*  Usually  E  is  expressed  as  a  function  of  q,p  =  ^4;  we  keep  q  rather  than  p  here. 
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Differentiation  by  time  gives 


•  _  dL  .  d  .  dL .  .  _  .  _ 

E  =  5 aj  +  qlt{ w’ "  ,q ~  “q  =  "*• Dii 

we  used  equation  (3.1)  in  the  last  step.  | 

Equation  (3.1)  has  the  same  invariance  property  as  the  conservative  Lagrange  equa¬ 
tions. 

Theorem  3.2. [12]  The  dissipative  Lagrangian  system  (S.l)  is  invariant  under  the 
change  of  variables  q  =  q(Q).  More  precisely ,  if  q{t)  satisfies  (S.l),  then  Q(t)  satisfies 
equations  of  the  same  form: 

jfq-Cq+V, j=0,  (3.1)’ 

Mere  C(Q,Q)  =  L(q(Q),q'(Q)Q)  and  D (<?,<?)  =  D(q(Q),q'(Q)Q). 

Proof:  A  simple  calculation  shows: 

—  £q  +  Vq  =  ^(Q)7 ^  ~  Lq  +  Dt'j 

I 
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4.  A  Rotating  Rigid  Body  with  a  Beam  Attachment 


Consider  the  spacecraft  depicted  in  Figure  4.1.  The  key  features  of  this  structure  are 
a  rigid  body  to  which  a  flexible  cantilevered  beam-like  appendage  of  length  l  is  attached. 


Figure  4.1:  A  rigid  body  with  cantilevered  beam  attachment. 

As  we  have  mentioned  in  section  2,  we  affix  the  moving  frame  to  the  rigid  body. 
More  specifically,  we  place  the  origin  of  the  frame  at  the  point  of  attachment  of  the  beam 
and  align  the  23  =  2-axis  along  the  undeflected  beam.  Viewing  this  cantilevered  beam  as 
essentially  a  one-dimensional  object,  we  describe  the  elastic  deformations  u(zs,t)  =  u(z,t) 
with  respect  to  the  coordinate  axes  (21,23,23)  =  (z,y, 2)  depicted  in  Figure  4.1.  More 
precisely,  u (z,t)  is  the  position  of  the  particle  whose  neutral  position  is  at  (0,0,2).  The 
decomposition  of  the  system  into  the  rigid  part  and  the  elastic  beam  corresponds  to  the 
decomposition  of  kinetic  energy  (2.1)  into  the  sum  of  rotational  and  translational  energies 
of  the  rigid  part  and  the  energy  of  the  beam.  (Note  that  u(-,<)  restricted  to  the  rigid 
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component  is  just  the  identity  mapping  for  all  t.)  We  have 


T  -  i uTIu  +  n v,yTYSle  +  |nu||y||2 

+U‘\\Y(nu  +  u,)  +  Hpdz 

where  S(ft)  =  u>,  c  is  the  center  of  mass  of  the  rigid  body  in  the  body  frame,  ms  is  the 
mass  of  the  rigid  body,  we  have  scaled  the  linear  mass  density  of  the  beam  to  be  one,  and 
the  inertia  tensor  with  respect  to  the  body  frame  is  given  by 

(  I*  Ixy  Ixx  \ 

I  —  Ixy  Iy  Iyt  j 

\/,z  Iyx  lx 

•where  Ix  is  the  moment  of  inertia  with  respect  to  the  z-axis,  etc.  (cf.  [1]). 

Since  the  beam  is  clamped  at  the  origin  of  the  (z,  y,  z)-coordinate  system  and  free  at 
its  other  end,  the  following  boundary  conditions  are  assumed: 

*(0,0  =  £s<0.0  =  £?«,«)  =  g(i.i)  =  0, ;  =  1,2, 

and  1*3(0, t)  =  =  0.  These  boundary  conditions  are  standard  in  the  theory  of 

ciamped-free  beams.  Here  Uj.vj.U’  are  the  deflections:  v  =  (111,1*2,2  +  1*3).  Note  that  1*3 

is  not  the  2-coordinate  of  u. 

The  equations  of  motion  for  our  rotating  satellite  are  obtained  by  incorporating  (2.5)- 
(2.7)  into  the  formalism  of  Lagrangian  mechanics,  as  discussed  in  the  previous  section. 
Thus  we  look  for  extremals  of  the  Lagrangian  L  =  T  —  V  with  kinetic  energy  T  given 
above  and  potential  V  given  by  the  strain  energy 

Viu)=ll‘ll“  W )!  +  "’M )’  +  «K)21  dz,  •  =  ■§;, 


17 


(4.2) 


where  only  quadratic  terms  were  retained  and  the  material  is  assumed  to  obey  a  linear 


Hooke’s  law.  Here  pi  (respectively  /z2)  gives  the  bending  elasticity  within  the  xz-plane  ( yz - 
plane  respectively),  and  /ij  is  the  Hooke’s  constant  giving  the  beam’s  stretching  elasticity. 
Unless  the  beam  is  abnormally  thick,  ^3  >>  pi,  /i2-  The  dissipation  function  is  given  by 


d = d(u,u)  =  \  j‘k,wr + MO’ + (4-3> 


where  il" ,  tij  can  be  thought  of  as  the  rates  of  change  of  appropriate  curvatures,  while 
U3  is  the  rate  of  change  of  the  contraction  coefficient  u'3.  The  ki’s  are  positive  constants 
reflecting  the  rates  of  energy  dissipation  due  to  deformation  of  material  in  the  beam. 


Remark  4.1.  It  might  seem  at  a  first  glance  that  the  first  two  terms  in  the  integral 
(4.2)  should  be  replaced  by  a  quadratic  form  (in  u'jtu!{) 


a(uj)2  +  2&u"  tx2  +  c(u2  )2; 


however,  by  properly  turning  the  body  coordinate  system  around  the  z-axis,  we  can  di¬ 
agonalize  this  form.  In  Figure  4.2,  the  r-axis  is  chosen  along  the  direction  in  which  the 
beam  bends  most  easily  (hi  <  p2) i  consequently,  the  beam  offers  the  stiffest  resistance  to 
bending  within  the  j/z-plane. 
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Remark  4.2.  We  must  point  out  that  the  beam  was  assumed  to  be  of  uniform 
crossection.  Expression  (4.2)  would  have  to  be  modified  to  include  the  beams  with  variable 
crossections  like  the  ones  show  in  Figure  4.3. 


Figure  4.3:  Beams  with  variable  crossection. 

The  modification  is  in  fact  quite  simple:  one  would  only  have  to  replace  the  first  terms  in 
(4.2)  by  a  quadratic  form  with  z- dependent  coefficients.  In  the  case  of  a  helical  beam  we 
would  take  the  matrix 

(:  ‘)  .*.«.>-*->(?  ;,)* 

where  R(z)  is  a  rotation  matrix  (  cof  a  8ina  j  ^-dependent  angle  a  =  const  •  z 
'  '  \  —  sin  a  cos  a  ) 

Remark  4.3.  To  incorporate  torsional  deformations  of  the  beam,  we  introduce  the 
torsion  angle  a(z;  t),  which  is  the  angle  formed  between  the  r-axis  and  the  projection  onto 
the  zy-plane  of  the  segment  rigidly  connected  to  the  beam  so  that  in  the  unperturbed 
position  of  the  beam  it  is  attached  at  (0,0,  z)T  and  is  parallel  to  the  s-axis.  Thus  the 
variable  a(z;  t)  describes  a  normal  bundle  of  the  beam.  Potential  energy  of  the  beam  is 
given  by 

V  =  i  /'(<  R-'„Rv",v"  >  +«(<ii)2  + 

2  Jo 

where  n  =  diag(fi  1,(12),  is  the  torsional  elasticity  coefficient,  and  v  =  (ui,U2). 


19 


A  slightly  more  subtle  remark  regarding  our  model  with  potential  energy  defined 
by  (4.2)  is  that  it  does  not  involve  potential  terms  incorporating  tensile  forces  into  the 
model.  Our  model  is  thus  somewhat  different  from  what  one  would  obtain  using  a  so- 
called  “geometrically  exact”  beam  theory.  It  is  perhaps  surprising  that  the  rotational 
equilibria  described  in  the  next  section  seem  to  differ  little  from  what  one  would  expect 
in  the  case  of  a  geometrically  exact  model.  We  refer  the  reader  to  the  recent  paper  [36] 
for  a  more  complete  discussion  of  the  effects  of  constitutive  restrictions  on  the  qualitative 
dynamics  of  rotating  systems. 

The  following  theorem  is  a  straightforward  consequence  of  the  results  presented  in  the 
previous  two  sections. 

Theorem  4.1.  Given  the  system  depicted  in  Figure  4-1  and  described  above  with 
kinetic  energy  (4-1),  potential  energy  (4-&)  and  dissipation  function  (4-S)>  the  equations  of 
motion  are  given  by 


Da  +  (mj,c  +  f  v.)  x  Y  1y  =  0 
Jo 


D2u  +  (idu  +  kdu  +  Y  1y  =  0 


ib(y  +  Yc)"+  [*YD2udz  =  0 
Jo 


where  the  quantities  in  these  equations  are  given  as  follows.  D(  )  =  £()+«x(),«(i)  = 
Iu(t)  +  J*  u  x  (Du)dz,  I  is  the  inertia  tensor  in  the  body  frame  defined  above,  m&  is  the 
mass  of  the  rigid  body  component,  Y  =  Y S-1(u;)  ( with  S(  ■)  is  as  defined  as  in  Lemma 
t.l),  and  d  is  the  differential  operator  defined  by  d  =  <?f),  n  =  diag{^i,p.2,n%) 

and  k  =  diag(ki,k2,k3). 


Equations  (4.4)-(4.6)  are  equivalent  to 


Iu>  +  u)  x  Iu}  +  /  u  x  [utt  +  «xu  +  2wxui  +  wx(«xu)  +  y  1y\dz 
Jo 


(4.7) 


+1114C  x  Y  *y  =  0 


Utt  +  u’xu  +  2wxuj  +  wx(wxu)  +  fidu  +  k&ut  +  Y~1y  =  0,  (4.8) 

m y  +  Y([  udz  +  mi,c)  =  c' +  c"t,  (4.9) 

Jo 

where  ms  denotes  the  mass  of  the  rigid  body  component,  and  m  denotes  the  total  mass 
of  the  body-beam  system.  Since  the  mass  density  of  the  beam  has  been  normalized  to  be 
1,  we  find  m  =  l  4-  m*.  (4.7)  can  be  further  rewritten,  using  (4.8),  as 

7u  +  ux/u—  f  u  x  [fidu  +  cdii\dz  +  m*c  x  T-1  y  =  0.  (4.7)' 

Jo 

One  easily  checks  that  the  total  angular  momentum  J(y+Yu)  x  £(y+Yu)dm  is  conserved, 
either  using  Noether’s  theorem  or  by  direct  computation. 


Remark  4.4.  If  the  center  of  mass  of  the  system  is  at  rest  with  respect  to  the  space 
frame,  we  may  assume  that  c'  =  c"  =  0. 


We  indicate  the  physical  meaning  of  the  terms  in  the  equations  (4.7)-  (4.9).  The 

sum  of  the  first  two  terms  in  (4.7)  is  interpreted  as  the  rate  of  change  of  the  rigid  body’s 
angular  momentum.  The  second  term  in  the  brackets  gives  the  inertial  force  on  the  beam 
due  to  the  body’s  angular  acceleration,  the  third  term  is  the  Coriolis  force,  the  fourth 
is  the  centrifugal  force  and  the  last  term  is  the  D’Alembert  force.  Thus  equation  (4.7) 
expresses  the  conervation  of  the  total  angular  momentum  in  space  expressed  in  the  body’s 
coordinate  system.  Equation  (4.8)  is  just  Newton’s  law  for  the  beam  expressed  in  the  non- 
inertial  body  frame  (the  D’Alembert  principle) — it  accounts  for  various  inertial  forces. 
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Equation  (4.9)  expresses  the  conservation  of  the  linear  momentum  of  the  whole  system.  It 
is  important  at  this  point  to  make  the  following  remark. 

Remark  4.5  There  is  an  apparent  paradox  associated  with  eq.  (4.7)':  one  might  expect 
to  be  able  to  express  the  integral  term  in  terms  of  u  and  its  derivatives  at  z  =  0,  since  the 
body  feels  the  beam  only  through  the  attachment  point.  As  it  turns  out,  this  expectation  is 
not  met.  In  formulating  our  continuum  mechanical  model  of  the  beam,  we  have  neglected 
certain  effects  such  as  torsional  deformations.  The  implicit  rigidity  in  our  model  leads  to 
this  nonlocal  character  in  the  equation. 
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5.  Asymptotic  Dynamics  of  a  Rotating  Elastic  Structure 


In  this  section  we  begin  an  analysis  of  the  asymptotic  behavior  of  the  body-beam 
structure  described  in  the  preceding  section.  For  finite  dimensional  dissipative  mechanical 
systems,  LaSalle’s  invariance  principle  [10]  can  be  used  to  show  that  states  asymptotically 
approach  a  minimal  invariant  subset  of  the  zero  set  of  the  dissipation  function  discussed 
in  Section  3.  In  Section  6,  it  is  shown  that  this  type  of  analysis  may  be  extended  to 
certain  infinite  dimensional  systems  with  features  in  common  with  our  body-beam  model 
described  by  equations  (4.4)-(4.6).  While  we  do  not  prove  that  all  solutions  to  (4.4)-(4.6) 
tend  to  the  zero  set  of  the  dissipation  function  D  as  defined  by  equation  (4.3),  we  do  offer 
a  more  or  less  complete  characterization  of  the  set  of  asymptotic  equilibrium  states  (by 
which  we  mean  the  set  of  solutions  to  (4.4)-(4.6)  for  which  the  equality  D  =  0  also  holds). 
We  prove,  in  particular,  that  in  asymptotic  steady  state  the  beam  displacement  function 
u(-  ,•)  does  not  depend  on  t.  Moreover,  it  is  shown  that  asymptotic  equilibrium  angular 
velocities  are  constant  vectors  parallel  to  the  principal  axes  of  the  steady  state  inertia 
tensor.  This  means  that  the  motion  of  the  system  is  a  pure  rotation  with  no  precession. 
This  is  the  content  of  the  following  theorem. 

Theorem  5.1:  Rrlcfir t  equilibrium  solutions  of  (4.7)-(j.9)  (by  which  we  mean  solu¬ 
tions  which  also  satisfy  D  =  0  with  D  as  defined  in  (4-3))  have  the  following  properties: 

(i)  there  i«  no  dependence  on  the  time  variable  t  in  the  beam  function:  u(z,t )  =  «oo(^); 

(ii)  the  angular  velocity  u;  is  a  constant  it's,; 

(Hi)  the  equilibrium  angular  momentum  is  a  constant,  ago  =  J ao<*>oo,  with  the  equilib¬ 
rium  inertia  tensor  of  the  combined  body-beam  system  given  by 

Joo  =  I+  f  uTuF.  -  uuTdz  -  m (ClCmE  -  cmc£)  (5.1) 

Jo 
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where  E  =  the  identity  matrix  and  Cm  =  ^(mjc  +  udz)  is  the  center  of  mass  of 
the  body -beam  system  (expressed  in  the  body  frame). 

(iv)  equilibrium  rotatuns  are  aligned  with  a  principal  axis  of  the  equilibrium  inertia 
tensor,  and  thus 

JooVao  =  Au;  OO  (5.2) 

Proof:  To  prove  (i),  note  that  for  equilibria,  J  ki(u")3  +  -f  ks(u'3)2dz  =  0. 

Hence  u"  =  u'f  =  tij  =  0,  and  the  result  follows  as  a  consequence  of  the  boundary 
conditions. 

Using  the  time-independence  of  u ,  we  show  that  corresponding  values  of  u>  are  con¬ 
stant.  Since  Ut  =  0,  equation  (4.8)  may  be  rewritten 


wxu  +  aix(wxu)  +  pdu  +  Y  1y  =  0. 


Differentiating  with  respect  to  z  at  z  —  0  and  using  the  boundary  condition  =  k  = 
0  J  at  t  =  0,  we  obtain 


iixk-(-wx(wxk)  +  fidux{0)  =  0 


Denoting  the  last  term  by  (ci,C2,C3)r,  rewrite  this  as 

(tii  —  U2W3  \  /  cj  \ 

wj  +wiwj  I  =  I  —Ci  J  .  (5-3) 

W1  +  /  \  C*  / 

Multiplying  the  first  two  components  by  and  wj  respectively  and  adding  we  obtain 


C2US1  —  CiU>2  =  0)2^2' 


(5.4) 


Prom  the  last  component  of  (5.3)  it  follows  that  the  right  hand  side  is  zero,  and  then  (5.4) 
together  with  the  last  component  of  (5.3)  implies  a>i  and  u>2  are  constant,  if  ci  and  C2  are 
not  both  zero. 
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To  establish  (ii)  in  the  case  that  c\  =  e2  =  0,  we  show  that  there  is  no  possible  choice 
of  body  coordinates  for  any  rotating  rigid  body  in  which  the  relationships  wy  =  u ;2o>3  and 
u>2  =  —  hold  as  in  (5.3)  unless  &U  o>j’s  are  constant.  We  proceed  by  noting  there  is  an 
orthogonal  change  of  basis  x  =  Bu ;  ( B  £  50(3))  such  that  the  angular  velocity  vector  x 
is  expressed  in  a  principal  axis  coordinate  system  with  components  satisfying  the  system 
of  differential  equations 

2i  —  aiXjXk 

where  j  =  k  =  cr2(x),  *  =  1)2,3,  and  <r  is  any  permutation  cn  the  three  symbols 

i  =  1,2,3.  It  may  be  shown  that  the  definition  of  the  coefficients  a*  in  terms  of  the 
principal  moments  of  inertia  implies  |oc, |  <  1,  and  ay  +  a2  +  <*3  +  ai 04013  =0.  The  form 
of  the  equations  u>j  =  u>2u>3,  <1>  =  —  ujyus  imposes  (polynomial)  conditions  on  the  a/s  and 
the  entries  of  the  orthogonal  matrix  B.  A  straightforward  analysis  shows  that  if  C  is  any 
orthogonal  matrix  and  if  the  transformed  vector  y  =  CTx  has  its  first  two  entries  obeying 
equations  of  the  form  yy  =  /Siyjt/j,  y2  =  PiUi  1/3 >  then  in  fact  C  must  be  a  permutation 
matrix  (a  nonsingular  matrix  in  which  each  column  has  a  single  1  and  all  other  entries 
equal  0).  The  /V s  therefore  must  be  equal,  modulo  a  permutation,  to  the  a^’s,  and  in 
particular  we  must  have  |/3j|  <  1,  which  is  not  consistent  with  the  form  of  the  equations 
for  wj,  u>2  unless  all  u;j’s  are  constant.  This  completes  the  proof  of  statement  (ii)  in  the 
theorem. 

(iii)  follows  from  a  direct  computation  using  (4.4)  and  (4.6).  It  may  also  be  obtained 
as  a  direct  consequence  of  the  parallel  axis  theorem. 

(iv)  follows  since  must  satisfy  u>oo  x  J 00^00  =  0.  | 

Corollary  5.1:  The  equilibrium  beam  function  u<x>(‘)  on<f  the  equilibrium  angular 
velocity  u>oo  related  by  equations  (5.1), (5.t)  together  with  the  fourth  order  system  of 
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ordinary  differential  equations 


ndu  =  -  cm),  (5.8) 

where  cm  is  the  center  of  mass  of  the  body-beam  system  in  the  body  frame,  u  =  (1*1,1*3,113  + 
2),  and  SCloo  =  Ugo,  and  where  the  boundary  conditions  are  as  prescribed  in  Section  4- 

Proof:  In  light  of  Theorem  5.1,  equation  (4.8)  may  be  rewritten  as 

u  =  -n^o  i*  -  y-1  y. 

Prom  equation  (4.9)  (in  which  c'  =  c"  =  0)  we  see  that  y(t)  =  —Y(t)Cm.  Moreover,  since 
y  satisfies  Y(t)  =  y(t)flOOJ  we  have  that  Y(t )  =  y0en~*.  Then  y  =  —  y(t)ft^0Cm  and 
Y~ly  —  —VllcCm,  proving  the  Corollary. 

Remark  5.1:  (5.8)  prescribes  a  nonlinear  boundary  value  problem  since  ftoo  depends 
on  i*oo(-)  through  equations  (5.1)  and  (5.2).  Some  idea  of  the  complexity  involved  in 
explicitly  determining  i*oo(-)  may  be  gleaned  from  our  solution  to  the  model  problem 
described  by  equations  (6.1)-(6.3)  in  the  following  section. 

Remark  5.2:  A  detailed  analysis  of  the  way  in  which  the  dissipative  dynamical 
system  (4.4)-(4.6)  evolves  toward  steady  state  is  beyond  the  scope  of  this  report.  Never¬ 
theless,  a  heuristic  description  which  makes  contact  with  the  classical  theory  of  rigid  bodies 
may  be  given  as  follows.  For  a  rotating  rigid  body  there  axe  two  well  known  conserved 
quantities:  the  kinetic  energy  E ,  and  the  magnitude  of  angular  momentum  \M\.  If  M 
represents  the  angular  momentum  vector  with  respect  to  the  body  frame,  then  the  energy 
E  =  1 MTI~1M  is  a  quadratic  form  in  M  with  constant  coefficients  and  the  conservation 
of  energy  confines  . 
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Figure  5.2:  Momentum  spheres  in  the  conservative 
and  dissipative  cases. 

M  to  an  ellipsoid,  resulting  in  figure  5.2a.  If  a  dissipative  elastic  appendage  is  present 
then  M,  while  still  confined  to  a  sphere,  will  move  to  smaller  and  smaller  energy  ellipsoids, 
as  indicated  in  figure  5.2b.  For  the  purpose  of  this  heuristic  discussion  we  ignore  the 
infinite- dimensionality  introduced  by  the  beam. 

It  should  be  noted  that  the  basins  of  the  two  sinks  Si  and  S2  in  figure  5.1b  are 
interlaced,  and  it  is  difficult  to  predict  which  sink  of  M  will  appear  if  M(0)  is  near  the 
maximum  energy  point  N.  This  somewhat  delicate  phenomenon  is  6aid  to  have  been 
overlooked  in  the  design  of  the  Explorer  II  mission.  The  satellite  oriented  itself  along  the 
pr  jper  axis,  but  in  the  direction  opposite  to  the  desired  one. 

It  must  also  be  noted  that  the  picture  in  figure  5.2  is  valid  for  moderate  \M\  only,  and 
it  undergoes  bifurcations  as  \M\  is  increased.  These  bifurcations  are  due  to  buckling  of 
the  beam  at  higher  angular  velocities.  Such  a  phenomenon  is  analyzed  in  complete  detail 
in  a  model  body-beam  problem  involving  a  1  d.o.f.  rotation  in  in  the  next  section  and  in 
a  model  problem  involving  a  mass-spring  attachment  to  a  spatially  rotating  body  in  [28]. 
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6.  Bifurcations,  Stability  and  Dissipation  in  a  Model  Problem 

In  this  section  we  illustrate  two  interesting  phenomena  in  a  model  rotating  body-beam 
system — bifurcations  and  the  existence  of  finitely  many  predetermined  angular  velocities 
which  are  the  asymptotic  limits  of  all  motions.  We  will  see  that  the  example  considered  be¬ 
low  is  an  infinite- dimensional  Morse  system.  It  is  similar  to  the  general  body-beam  system 
considered  in  the  preceding  sections.  We  also  point  out  an  unexpected  connection  between 
this  problem  and  that  of  studying  bound  states  of  a  nonlinear  Schrodinger  equation,  c.f. 

[151- 

Figure  6.1  shows  our  model  consisting  of  a  disk  with  a  beam  attached  to  its  center 
and  perpendicular  to  the  disk’s  plane. 


Figure  6.1:  A  model  problem  with  1  rotational  d.o.f. 


The  disk  is  constrained  to  rotate  around  the  z-axis  without  friction,  so  that  the  angular 
momentum  of  the  system  is  conserved.  The  beam  is  constrained  to  the  z-u-plane,  and  all 
the  deflections  are  parallel  to  the  u-axis.  Equations  of  motion  of  the  beam,  including 
internal  damping,  are 

putt  +  HV-xtxz  +  huxxxxt  =  u ;2u,  (6.1) 


pu7dz  +  I)  =  M, 


(6.2) 
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u(0,f)  =  u,(0,f)  =  u„(l,()  =  «,„(!,  i)  =  0, 


(6.3) 


where  p,  k  charactrize  elasticity  and  damping  respectively,  and  equation  (6.2)  expresses 
conservation  of  angular  momentum,  I  being  the  disk’s  moment  of  inertia.  From  now  on 
we  set  p  =  n  =  1  to  simplify  notation.  Steady-state  solutions  will  satisfy  the  ordinary 
differential  equation 

-utxxx  +u>2u  =  0.  (6.4) 

(6.4)-(6.3)  define  a  nonlinear  eigenvalue  problem  with  M  as  a  parameter,  the  nonlinearity 
lying  in  the  u-dependence  of  u>  in  (6.2). 

Remark  6.1.  The  stationary  problem  (6.4)-(6.3)  admits  a  revealing  variational  for¬ 
mulation:  its  solutions  are  the  critical  points  of  total  energy 


when  the  angular  momentum  M  is  fixed,  as  one  easily  checks.  This  has  a  simple  physical 
explanation:  As  the  beam  vibrates  energy  is  dissipated.  Ultimately  the  virbatioins  axe 
damped  out,  and  the  total  energy  of  the  system  is  extremized.  Actually,  for  some  excep¬ 
tional  (in  a  sense  which  we  make  precise  in  Theorem  6.1)  initial  conditions  the  limiting 
energy  value  will  be  not  minimal,  but  rather  critical.  (See  Theorem  6.1  and  Figure  6.4.) 
We  note  also  that 

E  =  -  f\u„t)2dz  (6.6) 

Jo 

Remark  6.2  on  Stability.  The  same  variational  formulation  suggests  a  stability 
criterion  for  the  dynamical  system  (6.1)-(6.3).  Namely,  consider  the  space  S  of  all  triples 
(u,u,u>)  with  angular  momentum  M,  i.e.  satisfying  (6.2).  If  the  stationary  solution  xo  = 
(u,0,u>)  minimizes  total  energy  (6.5)  on  S,  then  this  stationary  solution  is  stable,  even 
in  the  undamped  case.  This  means  that  for  all  initial  data  (u,ut,u;)t=o  in  S  sufficiently 
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close  to  the  minimal  (u,0,u>)  (in  the  energy  metric)  the  solution  will  stay  close  to  the 
minimizing  solution  xq  for  all  time.  This  is  quite  similar  to  the  standard  minimal  energy 
criterion  of  Lagrange;  here,  however,  the  energy  is  minimal  only  subject  to  the  angular 
momentum  being  constrained.  This  minimality  »  guaranteed  if  the  following  condition 
holds:  Ezx  +  A M„  >  0  on  the  tangent  space  to  M  =  const,  at  the  point  zo,  where  A  is 
the  Lagrange  multiplier:  Mx  =  A Ex . 

This  is  the  idea  behind  the  Energy-Casimir  method,  which  is  the  second  derivative 
test  in  the  presence  of  the  constraints  and  which  is  applicable  in  the  infinite-dimensional 
settings.  This  method  has  been  used  to  establish  stability  in  a  variety  of  problems  (see, 
e.g.,  the  work  and  references  in  Krishnaprasad  and  Marsden  [8]). 

A  complete  picture  of  the  behavior  of  our  model  system,  in  particular  of  the  global 
phase  portrait  and  its  bifurcations,  is  given  by  Theorems  6.1  and  6.2  below.  A  perhaps 
surprising  consequence  of  this  result  is  that  the  system  selects  a  discrete  set  of  angular 
velocities  for  its  stationary  rotations  independent  of  the  angular  momentum  M\  This 
phenomenon  is  observed  in  a  number  of  model  problems,  but  it  is  by  no  means  universal. 
See  [36]  on  this  point.  Another  aspect  of  this  result  is  that  our  model  is  an  infinite¬ 
dimensional  Morse  system. 

To  state  our  first  result  we  combine  (6.1)  and  (6.2)  into  a  single  nonlinear  equation 

««  +  «**«  +  kuXSMMt  =  u>3(u)u 

where  u;(u)  =  M/(\\u\\2  +  /),  lju||2  =  J*  u3dz.  This  can  be  rewritten  as  a  system 

Theorem  6.1.  Qualitative  behavior  of  the  system.  Assume  that  k  >  0,  i.e.  the 
damping  is  present. 
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There  exists  an  infinite  sequence  0  <  u\  <  u>2  <  . . .  — ♦  oo  of  preferred  angular 

velocities  associated  with  the  problem  (6.7)  with  boundary  conditions  (6.S)  such  that  if 

the  angular  momentum  M  lies  in  the  interval  (M*,  Mfc+i)  =  (Jo;*,  Iwi,+i),  the  problem 

(6.7), (6. S)  has  k  distinct  nontrivial  stationary  solutions  (u,t>)  =  (11,(2),  0),  1  <  *  <  k,  and 

a  trivial  solution  (u,t>)  =  (0,0),  with  angular  velocities  «*>(«<)  =  — pf - =  <*>«,  1  <  i  <  k 

'Jo  * 

and  u;(0)  =  ™  These  solutions  are  depicted  in  Figure  6.t.  The  solution  with  the  smallest 
angular  velocity  is  dynamically  stable,  i.e.  any  solution  of  (6.7), (6. S)  with  nearby  initial 
conditions  remains  close  to  it  for  all  time,  in  the  LP-norm  given  by 

\\(u,u)\\  =  \\(u,v)\\=  [  (u 2  +  v7)dz. 

Jo 

Stationary  solutions  with  higher  angular  velocities  are  unstable,  and  furthermore,  i-th  so¬ 
lution  (in  the  order  of  increasing  angular  velocity)  (u,  v)  =  («i,0)  has  (»  —  1)- dimensional 
unstable  manifold;  this  dimension  coincides  with  the  number  of  zeroes  of  u,.  In  particular, 
if  M  <  Mi  =  Iuii ,  only  one  mode  u  =  0  is  present  and  is  stable.  As  M  crosses  the  bifur¬ 
cation  values  M{,  the  system  undergoes  a  series  of  pitchfork  bifurcations  shown  in  Figure 
6.3. 


u 


u 


2 


U0  =  0 


-u 


2 


-u 


I 


Figure  6.2:  Solutions  to  (6.1)-(6.3) 

There  have  been  several  recent  studies  of  the  evolution  of  solutions  to  infinite  dimen¬ 
sional  systems  toward  equilibria  (see  Dafermos  [4]  and  references  therein).  Theorem  6.2 


31 


below  shows  that  under  mild  (physically  reasonable)  assumptions  solutions  of  our  systems 
approach  the  equilibrium  states  described  above. 


Theorem  6.2  Fix  M  >  0  and  let  k  >  0,  May  k  =  1.  .Any  solution  ti>(<)  of  (6.7)  with 
xto(-z),  v0(^)  €  J5T(12)(0,1)  tends  in  the  L2-norm  as  t  — »  oo  to  one  of  the  stationary  solutions 
from  Theorem  6.1  (See  Figure  6.4). 


Figure  6.3  s  denotes  stable  equilibria,  and  Uj,  Uj, . . .  denote  unstable  equilibria 
with  1,2,...  dimensional  unstable  manifolds.  The  numbers  and 
I?  - norm  of  equilibria  depend  on  the  parameter  M. 


Figure  6.4.  Global  phase  portrait  of  the  system  (6.7), (6. S).  Indices  show 
the  dimension  of  the  unstable  manifolds,  which  are  depicted  by  arrows. 
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Proof:  of  Theorem  6.1. 


1.  Existence  of  k  +  1  solutions  for  Af*  <  M  <  Mk+i  •  Stationary  solutions  satisfy  the 
nonlinear  boundary  value  problem  (6.3)-(6.4).  There  exists  a  sequence  0  <  u>\  <  u>2  <  . . .  of 
angular  velocities  for  which  this  problem  has  nontrivial  solutions  (the  eigenfunctions  of 
with  boundary  conditions  (6.3))  which  we  denote  by  Cje<,  c,-  arbitrary  and  IM*)IU»  =  i- 
Define  Mj  =  Iuij*  Freedom  of  choice  of  c;  is  used  to  satisfy  the  angular  momentum 
constraint  (6.2),  which  becomes 

u>i(c2i  + 1)  -  M. 


This  can  be  solved  for  Cj  if  and  only  if  M  >  Iu){  =  Af,-;  setting  ti,  =  c,ej  proves  the 
existence  of  k  solutions  if  M  >  M*.  One  additional  solution  is  obtained  by  setting  Ufc+i  = 


0,  w*+i  = 

2.  Stability  of  stationary  modes.  We  restrict  ourselves  to  the  undamped  ( k  —  0)  case, 
as  the  dissipation  does  not  affect  stability  of  stationary  modes.  The  governing  equations 
can  be  written  as 


u-tt  —  ~nXZzz  "b 


(6.8) 


The  dimension  of  the  unstable  manifold  of  the  mode  u*(.z)  is  given  by  the  number  of 
positive  eigenvalues  of  the  linearization  of  the  operator  >l(u)  at  u  =  ttj.  This  linearization 
is  given  by 


A'(u)v  =  —d*v  + 


4 M2 

(ll«||J  +  I)s 


{utv)u, 


*  Mj  is  the  total  angular  momentum  when  the  beam  is  undefiected,  since  the  straight 
beam’s  (u  =  0)  contribution  to  the  angular  momentum  is  zero. 
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where  d  =  £  and  (u,t>)  =  J*  u(z)v(z)dz. 
At  u  =  Ui  we  have 


Bv  =  A'(ui)v  =  —  d*v  +  u>*v  —  M  1(u<,v)i*i 


The  last  term  in  this  expression  is  a  scalar  multiple  of  the  orthogonal  projection  of  v  onto 
Ui.  Comparing  B  with  the  operator  C  given  by 


Cv  =  —d*v  +  ufjv, 


we  see  that  the  spectrum  of  B  can  be  obtained  from  that  of  C  by  replacing  the  eigenvalue 
0  of  C  corresponding  to  the  eigenfunction  u,-  by  — 4u>*M_1  ||u»||2.  Thus  the  eigenvalues  of 
B  are 


2  2 
>ui  “  "i+iv 


of  which  precisely  t  —  1  are  positive.  This  shows  that  the  dimension  of  the  unstable  manifold 
at  the  stationary  solution  Ui  is  (*  —  1). 


Proof  of  Theorem  6.2  consists  of  two  steps: 

(1)  Showing  that  for  any  e  >  0  the  solution  w{t)  enters  and  stays  in  an  e-neighborhood 

(in  L 2  norm)  of  the  zero  set  of  the  dissipation  function.  This  set  is  {(u,t>)  :  v\tdx  = 

0}  =  {(u,v)  :  (u(®),0),  u  €  see  (6.6). 

(2)  Showing  that  if  w(t)  stays  in  an  e-neighborhood  of  the  zero  set  of  {t;  =  0},  then 
w  is  e-close  to  one  of  the  stationary  points  of  the  flow,  figure  6.5. 
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Figure  6.5  Solutions  which  remain  e-close  to  {t>  =  0}  must  be 
e-close  to  a  stationary  point. 

Proof  of  Step  1.  Pick  e  >  0,  and  consider  two  neighborhoods  of  {v  =  0}:  Afe  = 
{(u,  v)  :  HU*  <  e}  and  Jff/2-  Our  solution  tu(t)  enters  Jft/2  for  some  t  >  0  with  necessity: 
otherwise  for  all  t  large  enough  we  would  have 

E  =  -  [  u\ttdz  =  -[  vltdz<-\0(  v2dz  <  A0(J)2,  (6.9) 

Jo  Jo  Jo  * 

resulting  in  E(t)  — >  —  oo,  which  is  a  contradiction.  Here  X0  is  the  smallest  eigenvalue  of 
d4  =  (d*/dz4)  with  boundary  conditions  (6.3). 

Thus  to  prove  that  w(t )  stays  in  Jft  forever  after  some  T  >  0  it  remains  only  to 
exclude  the  possibility  of  infinitely  many  trips  between  Jft/2  and  the  exterior  of  Jft,  see 
figure  6.5.  We  will  do  so  by  showing  that  each  trip  results  in  a  loss  of  energy  at  least 
A E  >  0  depending  only  on  e  and  w(0)  at  each  crossing,  so  that  w(t)  can  afford  only  finite 
number  of  trips. 

Let  <i,<2  be  two  consecutive  crossing  times  of  the  boundaries  of  Afe  and  Jff/2-  We 

have 

\E(t,)  -  =  \  f‘’  (6-10) 

where  we  have  used  the  fact  that  J*  v2dz  >  (|)J  and  the  Poincare  type  estimate  as  in 
(6.9).  It  remains  only  to  provide  a  lower  bound  on  the  trip  time  |<a  —  fj  |;  it  is  provided  by 
the  upper  bound  on  the  velocity  w(t)  which  is  the  result  of  the  smoothness  assumption. 

Lemma  6.1.  Any  solution  w(t)  =  (u(z,f),v(r,f))  of  (6.7)  with  u 0(z),v0(z)  6  1) 

has  velocity  bounded  in  the  L 2  norm:  i.e.  there  exists  C  =  C( u0,v0)  >  0  such  that  for  all 
t  >  0 

||F(ti;(f))||2  =  HU,  4-  ||  -  d*u  -  9*v  +  ^(u^lU,  <  C2. 
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Proof  of  this  lemma  is  given  in  the  Appendix. 

By  the  choice  of  <i,<2  and  using  Lemma  6.1,  we  obtain 


|  <  M<»)  - -»(<.)!!  =  II J” f(.»W)<k||  <  |f,  -  «,| 


8up||F||  <  C\t7  -tj|, 

t>  0 


implying 


1*2  “ 


which  together  with  (6.10)  proves  that  w(t)  stays  in  Aft  for  all  t  >  T(e). 


Proof  of  step  t.  We  show  now  that  if  w(t)  stays  in  Aft  for  all  t  >  T(e)  then  it  must 
tend  to  an  equilibrium  point  of  the  flow  (6.7)  -more  precisely,  we  will  show  that  for  every 
S  >  0  there  exists  an  e  >  0  such  that  if  ||v||  <  e  for  all  t  >  T(e)  then  for  some  equilibrium 
solution  we  =  (u*,0)  we  have  ||ue  —  tt||x,a  <  6  for  all  t  >  T(e).  Our  strategy  is  to  show 
that  if  a  solution  is  u(<)  not  close  to  an  equilibrium,  then  its  velocity  v(t)  must  grow  thus 
taking  it  outside  the  neighborhood  Aft,  leading  to  a  contradiction. 

Assume  the  contrary:  There  is  some  6  >  0  such  that  for  every  e  >  0  there  is  a 
ta  =  tn(e)  >  T(e)  such  that  although 


||u(t)|!  <  e  for  all  t  >  T(e) 


(6.11) 


we  have 

IK<„)  -  «.||  >  «.  (6.12) 

for  every  equilibrium  ue.  The  idea  is  to  show  that  being  far  from  any  equilibrium  causes 
an  increase  in  velocity  thus  causing  ||v||  >  e  (a  contradiction).  The  details  are  as  follows. 
(6.12)  implies  for  f{u)  =  —  d*u  +  u;2(u)u: 

||/(tt(/0))||>A||«(<0)-Tie||>A«,  (6.13) 
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for  some  A  >  0  independent  of  6 ,  if  we  choose  6  small  enough  (which  we  d.S). 

We  have:  (cf.  (6.7)) 

ll»(<)  -  »(<.)IU-  >  II  jf*  /(<*W)<k||  -  II  jf  fi*»(*)i»||.  (6.14) 

From  the  proof  of  Lemma  6.1  it  is  clear  that  there  exists  C  >  0  such  that  ||i*  <  Ce 
for  t  large  enough;  without  loss  of  generality  we  assume  that  this  already  holds  for  t  >  t0. 
We  have  from  (6.14): 

M*)\\  >  II  T  fM)M  -  <1  +  C)(t  -  t0)  for  all  t  >  t0  (6.15) 

Jt. 

Furthermore,  if  t'  ,t"  are  sufficiently  large,  and  if 

Hi")  -  «(t')||  <  v? 


holds,  then 

||/(«(i"))  -  /(«(!'))||  <  l||/(«(i'))|| 

This  can  be  seen  from  the  proof  of  Lemma  6.1  which  shows  that  all  the  harmonics  of  suffi¬ 
ciently  high  order  decay  exponentially,  with  w(t)  tending  to  a  finite-dimensional  subspace 
of  X 2  x  I-2.  Applying  this  remark  to  (6.15)  with  t1  =  t0,t"  =  t,  we  obtain  the  final  lower 
bound  on  v,  using  (6.13): 

IM«)II  >  j(<  -  «.)ll/(»(«.))ll  -  <(1  +  C)(<  -  *.)  >  (6-16) 

(i-l.)|iM-e(J  +  C)], 

as  long  as 


||u(<)  -  u(<o)||  <  y/e. 


Using  this  estimate  on  the  velocity  v(t),  we  will  show  that  w(t)  leaves  the  box  ||t>||  < 
e,  jju  —  ti(t0)||  <  \/e  through  the  “horizontal”  boundary,  i.e.  ||v(<*)||  =  e  for  some  first  exit 
time  t*  >  ta,  see  figure  6.5.  This  will  result  in  the  desired  contradiction. 
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To  see  that  ||v||  =  e  is  reached  first,  we  note  that  otherwise 


||u(<*)-u(*o)||  =  v'e, 


and 

e(<*  -  to)  >  ||  /*  *(0*||  =  || *(«•)  -  u(io)||  = 

Jt. 

i.e.  the  time  required  to  reach  the  vertical  boundary  is  large: 


tm  -t0  > 


V*' 


t* 

Using  this  in  (6.16) — which  is  valid  for  ta  <t  <* - we  obtain 


-A5-e(l  +  C) 


>  e, 


which  is  a  desired  contradiction,  if  we  pick  e  sufficiently  small.  {I 

This  shows  that  given  any  8  >  r  there  exists  an  e  >  0  such  that  if  w  stays  in  .A/’t,  then 
u(z,<)  is  less  than  8  (in  Z2[0,1])  from  an  equilibrium  solution  tie  of  (6.8).  Since  the  first 
step  in  this  proof  showed  that  for  e- neighborhood  of  the  zero  set  {v  =  0}  of  the  dissipation 
function,  this  proves  Theorem  6.2. 
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|0,!»o||  =  X>J(0)Ai  <  c 


imply 

M°)l>  |4*(0)|  <  cA*s. 


We  will  show  that  for  all  t  >  0  the  estimates 


M<)li  \bk(t)\  <  c\k  2 


hold,  for  k  large  enough,  thus  implying  (A. 2): 

£  ^ lal  -  ^2  4  =  g2  ^k  2  <  °° 


since  r^i  k4. 

To  show  (A. 3)  we  rewrite  (A.l)  as 


z  =  Az  ■+•  Rz 


where 


and 


Introducing  the  Lyapunov  matrix 

B  =  f  eA  *eAtds  satisfying  ATB  +  BA  =  —I, 
Jo 


we  obtain 

=  ~(z,z)  +  ((RTB  +  BR)z,z)  < 


-(1 -||«T£ +  £*!!)(*,  *)■ 


(A.Z) 
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An  explicit  computation  gives 


’I  +  J-  J_ 

1  ^  2A  2A 


B  = 


J_  J_  ,  1 

2A  2A  'r  2A* 


and  thus  || RTB  +  BR\\  <  since  u  <  *f- 

We  have 

<£  .  _  .  .  “const."..  .  1.  . 

—  (Bz,z)  <  -(1 - - - )(*,*)  <  --(*,*), 

the  latter  inequality  holding  for  all  A  =  A*  large  enough. 

Now,  the  minimal  eigenvalue  p  of  the  positive  definite  matrix  B  is  estimated  from 
below  as 


P>2A 


(for  A  large  enough). 
We  have  for  all  t  >  0  and  for  k  large  enough 


2A 


or 


(z(t),z(f))  <  (£z(f),z(t))  <  ( BzoyZo )  <  (*o,.zo), 

al(t)+bl(t)<2Xk(al(0)  +  bl(0))< 

<  2A*caA*8  =  2caA*5. 


This  implies  (A.3),  q.e.d. 
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EQUILIBRIUM  MECHANICS  OF  ROTATING  SYSTEMS 


John  Baillieul 
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Boston.  MA  02215 

ABSTRACT 

Equilibrium  rotation#  of  a  planar  kinematic  chain  undergoing  a  aingle  degree  of  freedom  (out  of 
plane)  rotation  in  a  gravitational  field  are  studied.  Sharp  bounds  on  the  number  of  dynamic  equilibria 
are  given,  and  the  dependence  of  this  number  on  physical  parameters  is  discussed  for  the  special  case 
of  two-link  chains  and  chains  with  identical  links.  The  theory  indicates  that  there  will  be  a  richer  set 
of  dynamically  stable  equilibria  for  carefully  designed  chains  with  unequal  links.  This  is  born  out  by 
laboratory  experiments. 


1.  INTRODUCTION 

As  a  prelude  to  a  comprehensive  study  of  the  dynamic  stability  of  rotating  mechanical 
systems,  it  is  useful  to  study  the  equilibrium  dynamics  of  various  classes  of  prototypes. 
Indeed,  a  number  of  participants  at  this  conference  have  carried  out  such  studies  involving 
rotating  mechanisms  featuring  various  different  types  of  elastic  and  articulated  appendages.  In 
our  own  recent  work  (Baillieul  and  Levi.  1987),  the  rotational  dynamics  of  a  simple  structure 
consisting  of  an  elastic  rod  attached  to  a  rigid  body  were  studied.  The  transient  dynamics  for 
structures  of  this  type  involve  vibrations  of  the  attachment  which  are  subject  to  inertial  forcing 
and  viscoelastic  damping.  We  have  shown  that  in  the  absence  of  external  forcing,  the  damping 
will  have  a  dominant  influence  on  the  long  term  dynamics,  and  ultimately  all  movement  of  the 
appendage  with  respect  to  the  rigid  body  coordinate  frame  will  disappear.  Moreover,  we  have 
shown  that  in  steady  state,  rotations  always  tend  to  a  constant  velocity,  and  for  the  three 
rotational  degree  of  freedom  problem,  this  angular  velocity  vector  is  aligned  with  a  principal 
axis  of  the  steady  state  inertia  tensor  of  the  body  beam  system.  While  we  do  not  yet  have 
a  complete  understanding  of  the  nonlinear  relationship  between  steady  state  angular  velocity 
and  steady  state  elastic  deformation  of  the  beam  appendage,  it  has  been  possible  to  carry 
out  a  detailed  analysis  of  the  qualitative  dynamics  of  a  closely  related  model  having  only 
one  rotational  degree  of  freedom.  In  this  system,  the  rotation  is  restricted  to  be  about  the 
neutral  axis  of  the  beam.  All  steady  state  solutions  and  corresponding  stability  characteristics 
have  been  determined,  and  an  explicit  description  of  the  system's  dependence  on  the  angular 
momentum  has  been  developed.  It  has  been  shown  that  as  angular  momentum  is  increased, 
new  equilibria  appear  as  a  result  of  a  sequence  of  pitchfork  bifurcations.  No  matter  how  rich 
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the  equilibrium  set  becomes,  however,  there  is  at  most  one  stable  equilibrium  configuration  if 
we  assume  the  beam  is  uniform.  (Below  a  certain  threshold  angular  momentum,  the  neutral 
or  undeflected  state  of  the  beam  is  stable,  while  above  this  threshold  the  undeflected  state 
ceases  to  be  stable,  and  the  first  mode  shape  becomes  the  stable  configuration  for  all  higher 
values  of  angular  momentum.) 

The  bifurcations  displayed  in  the  rotating  beam  model  are  similar  to  phenomena  described 
in  classical  works  on  whirling  shafts  (see  e.g.  Dick.  1948).  In  order  to  obtain  a  more  detailed 
understanding  of  this  type  of  buckling  in  the  context  of  the  modern  qualitative  theory  of  non¬ 
linear  dynamics,  we  have  turned  our  investigations  to  the  analysis  of  rotating  planar  kinematic 
chains.  In  the  next  section,  we  provide  a  precise  statement  of  the  problem  under  investigation 
and  summarize  recent  results  (detailed  in  Bailiieul.  1987)  regarding  the  numbers  of  equilibria. 
In  Section  3.  we  present  new  results  which  suggest  that  chains  which  are  uniform  (i.e.  have 
nearly  identical  links)  will  have  a  less  rich  set  of  interesting  dynamically  stable  equilibrium 
configurations  than  chains  in  which  the  links  have  carefully  prescribed  unequal  masses  and 
lengths. 

2.  EQUILIBRIA  IN  A  ROTATING  CHAIN 

In  this  section  we  consider  the  rotational  dynamics  of  the  rotating  simple  kinematic  chain 
depicted  in  Figure  2.1.  Here  a  certain  number,  n.  of  links  of  various  lengths  are  connected  to 
form  a  single-strand  planar  kinematic  chain.  We  suppose  this  chain  is  suspended  in  such  a  way 
that  the  force  of  gravity  tends  to  extend  the  chain  to  its  maximum  total  length.  We  further 
assume  there  is  (360°)  free  rotation  (with  no  actuation  or  friction)  about  all  joints  in  the  planar 
chain,  but  that  torque  may  be  applied  about  the  vertical  axis  (which  passes  through  all  links 
and  joint  axes  when  the  linkage  is  in  the  neutral  (fully  extended)  configuration).  We  suppose 
the  links  are  massless,  but  each  joint  in  the  planar  chain  has  mass  m;_i  (i  =  2, . . .  ,n).  and 
there  is  a  mass  m„  located  at  the  tip  of  the  final  link.  The  angles  V’t  depicted  in  Figure  2.1 
measure  deviation  of  each  link  from  its  neutral  vertical  position.  We  also  let  o  denote  the 
angle  of  rotation  of  the  mechanism  about  the  vertical  axis.  (V»o  is  measured  with  respect 
to  an  arbitrarily  chosen  reference.)  Depending  on  the  mass  distribution  in  the  links  of  an 
actual  physical  chain,  the  mathematical  idealization  we  have  chosen  to  study  may  or  may  not 
provide  a  model  which  will  predict  detailed  dynamical  behavior.  The  advantage  of  working 
in  this  simplified  setting,  however,  is  that  we  are  able  to  display  a  wide  range  of  qualitative 
features  and  to  understand  the  way  in  which  these  depend  on  parameters. 
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Figure  2.1:  An  n- degree  of  freedom  planar  kinematic  chain  in  a  gravitational  field  rotated 
about  its  neutral  axis. 

Following  a  geometric  formalism  along  lines  similar  to  those  detailed  in  Baillieul  and  Levi 
(1987).  it  is  strightforward  to  write  down  the  dynamical  equations  of  this  rotating  chain.  At 
present,  however,  we  are  simply  interested  in  the  special  case  where  =  u  (a  constant)  and 
ipj  —  0  for  j  =  1, . . .  ,n.  Such  equilibrium  solutions  to  the  dynamical  equations  are  also  given 
as  the  zeros  of  the  gradient  system 


where 

L  = 

=  sinV>i)2  +  mj(£i  sinV»i  + 1-2  sin^a)2  + 

sinV’i  +  ...  +  lnBinr(>n)2}v2 
—  [mi£icosV»i  +  fua(/i  cos^i  +^acos^j)+ 

•••  +mn(£icosil>i  +  ...  +  lncoBiJ>n)}g 

with  g  designating  acceleration  due  to  gravity.  This  system  of  equations  may  be  explicitly 
written  out 
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[(mi  +  .  .  .  4-  mn)liSi  +  (m2  +  .  .  •  +  Tnn)l282  +  •  ■  •  +  mn£nSn]ciU;2 
-(mi  +  ...  +  mn)aig  =  0 


[(m2  +  . .  •  4-  mn)^i  Ji  +  (m2  4- ...  4-  mn)t2a2  4-  . . .  4  mn£n£n]c2t*ja 


— (m2  4- ...  4-  mn)»2g  =  0 


Mi»i  +  rnnl2a2  +  . . .  +  mnlnan]cnu2  -  mns„g  =  0.  (2.1) 

where  as  usual  =  sink’s  and  c*  =  cos  V»». 

Frictional  forces  will  oppose  motions  in  the  joints  of  these  kinematic  chains,  and  the 
stability  of  each  equilibrium  will  be  characterized  by  the  signs  of  the  eigenvalues  of  the  Hessian 
|^y.  Before  discussing  the  stability  characteristics  of  various  equilibria,  we  establish  upper 
and  lower  bounds  on  the  numbers  of  solutions  to  (2.1).  The  exact  number  of  solutions  in  any 
particular  case  will  be  seen  to  depend  on  the  values  of  the  parameters  m*,  and  u>.  It  will 
be  shown  in  the  next  section  that  for  the  case  n  =  2  the  bounds  we  have  obtained  are  tight. 


Theorem  2.1:  For  small  values  of  ju;j,  there  are  precisely  2n  distinct  solutions  to  the 
system  (S.l).  These  solutions  are  given  explicitly  by  {(^i,...,^n)  :  =  0  or  7r}.  The 

index  ( =  number  of  negative  eigenvalues  of  the  Hessian  ))  °f  anV  *uch  solution  is  the 
number  of  vector  entries,  ipj,  which  are  equal  to  ir. 

This  theorem  is  proved  in  Baillieul  (1987).  Physical  intuition  suggests  that  there  will 
be  more  than  2n  equilibrium  solutions  as  the  parameter  u>  (or  equivalently  /)  is  increased 
in  (2.1).  To  obtain  a  sharper  picture  of  the  way  in  which  bifurcations  occur,  we  simplify 
the  mathematics  by  means  of  some  assumptions  and  a  change  of  variables.  First,  note  that 
for  the  purposes  of  establishing  a  qualitative  theory  of  solutions  to  (2.1).  there  is  no  loss 
of  generality  in  normalizing  g  to  be  1  and  only  allowing  u>  and  the  l<,  m,  parameters  to 
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vary.  Then  letting  Xi  =  cos  ipi.  yi  =  sin  (2.1)  is  transformed  into  a  system  of  algebraic 
(polynomial)  equations: 


[(mi  +  •  •  •  +  Tnn)£iyi  +  (mj  +  ■  •  •  +  rnn)£22/2+ 


•  •  •  +  m„inyn]x iu2  -  (mi  + - |-  m„)yi  =  0 


[£l3/l  +  £21/2  + - 1-  ^nVn]*n«J  ~  1/n  =  0 

xl  +  y\  ~  1  =  0 


xn  +  y»  ~  i  =  0  (2.2) 

We  shall  call  the  system  (2.2)  the  algebraic  equilibrium  equations  for  the  rotating  planar 
chain.  The  precise  connection  between  solutions  to  (2.1)  and  (2.2)  is  given  in  the  following 
easily  proved  lemma. 

Lemma  2.1:  Real  solutions  to  (t.l)  define  the  solutions  to  (t.l). 

While  there  are  2n  equations  in  2n  unknowns  in  (2.2) — a  seeming  increase  in  complexity 
over  (2.1) — we  may  appeal  to  classical  algebraic  geometry  and  intersection  theory  in  deter¬ 
mining  the  number  of  solutions  to  (2.2).  Indeed,  since  there  are  2 n  quadratic  equations. 
Bezout's  theorem  suggests  there  will  be  as  many  as  2Jn  =  4n  (possibly  complex)  solutions. 
The  precise  number  is  given  by  the  following  theorem,  whose  proof  is  given  in  Baillieul  (1987) 
without  reference  to  the  classical  theorem  of  Bezout. 

Theorem  2.2:  For  every  choice  of  angular  velocity  u>  and  mass  and  link-length  param¬ 
eters  mi, . . .  ,m„,  £i,.  ..,£„,  there  are  4n  (some  possibly  complex)  solutions  to  the  system 


The  proof  here  depends  on  the  following  lemma  (also  proved  in  Baillieul.  1987)  which  is 
of  interest  in  its  own  right. 

Lemma  2.2:  As  u>  — *  oo,  solutions  to  (t.t)  remain  finite. 
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Because  there  is  a  one-to-one  correspondence  between  solutions  to  (2.1)  and  real  solu¬ 
tions  to  the  corresponding  algebraic  equilibrium  equations  (2.2),  Theorems  2.1  and  2.2  taken 
together  imply  the  following  result. 

Corollary  S.l:  Given  any  values  for  the  parameters  there 

are  at  least  2n  but  no  more  than  4n  solutions  to  the  equilibrium  equations  (S.l). 


3.  THE  ROTATING  DOUBLE  PENDULUM  AND 
THE  CASE  OF  UNIFORM  CHAINS 

The  analysis  of  the  preceding  section  shows  that  for  a  two  link  chain  there  are  between  4 
and  16  equilibrium  states.  To  facilitate  a  detailed  analysis  of  the  dynamics  of  this  mechanism, 
we  reformulate  equations  (2.1)  by  letting  g  =  1,  t  =  t\,  l2  =  rt,  (0  <  r  <  oo),  a  = 
m2 /(mi  +m2).  Under  these  assumptions  (2.1)  may  be  rewritten 

(^sin^i  +  art  sin  if  2)  cos  ifiu>2  —  sin  if\  =0 

(^sin^i  +  rising)  cos  —  sin  ^2  =  0-  (3-1) 

The  following  theorem  gives  a  substantial  answer  to  the  question  of  parametric  dependence 
in  the  case  n  =  2. 

Theorem  S.l:  Let  a  and  r  be  fixed  in  (S.l),  and  suppose  r  >  1. 

(i)  As  uj2  increases  from  0,  four  pitchfork  bifurcations  occur  involving  the  fixed  equi¬ 
librium  points  (V’l.V’a)  =  (0, 0),  (rr, 0),  (0,7r).  For  sufficiently  large  values  of  u>2,  these 
bifurcated  equilibria  account  for  twelve  of  the  possible  solutions  to  (S.l). 

(ii)  If  a  bifurcation  other  than  accounted  for  in  (i)  occurs,  it  involves  the  appeamce  of 
a  symmetric  pair  of  saddle  nodes. 

(Hi)  If  ar  >  1,  then  for  sufficiently  large  values  of  w2  there  are  only  twelve  solutions 
to  (S.l).  Two  are  stable  equilibria  (the  index  of  the  Hessian  of  L  is  0),  six  are  index  1 
critical  points  of  L,  and  four  are  index  t  critical  points.  If  ar  <  1,  then  for  all  sufficiently 
large  values  of  w2  there  are  sixteen  solutions  (counting  multiplicities)  to  (S.l).  Four  are 
stable,  eight  have  index  1,  and  four  have  index  i. 
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Proof:  (i)  It  is  an  elementary  calculation  to  show  that  the  Hessian  of  L  is  represented 
by  the  2x2  matrix 

/  ,  x  / ^2[(*1  +  ar52)3i  -  c\ ]u>2  +  ic  1  -ar£2dc2u2 

'  1  2'  ^  —arl2c\c2u)2  ar£2[(si  +  rs2)»2  —  rcl]u2  +  ar£c2 

Bifurcation  of  solutions  to  (3.1)  occur  when  this  matrix  becomes  singular.  For  the  equilibrium 
points  (0,0),  (7r,0),  and  (0,7r),  the  Hessian  has  repspective  singular  loci  satisfying 

r(l  —  a)  u>*  —  (1  +  +1  =  0 


r(l  —  a)u;4  +  (r  —  l)u>2  —1  =  0 
r(l  —  a)w4  +  (1  —  r)u )2  —1  =  0. 

(The  matrix  is  never  singular  for  the  unstable  equilibrium  (*,*)■)  Noting  that  only  real 
solutions  of  these  equations  are  of  interest,  we  find  they  may  be  rewritten  to  provide  formulae 
expressing  u>2  as  a  function  of  r  and  a: 


(a)  u,2 


r  +  1  ~  y/(r  +  l)2  -  4r(l  -  a) 
2r(l  -  a) 


(6)  a,2 

(c)  u,2 

( d )  w2 


1  -  r  +  y/(l  -  r)2  +  4r(l  -~^) 
2r(l  —  a) 

r  -  1  +  y/(r  -  l)2  +  4r(l  -  a) 
2r(l  —  a) 

r  +  1  +  y/jr  +  l)2  ~  4r(l  -  a) 
2r(l  —  a) 


These  expressions  for  a.’2  are  listed  in  order  of  increasing  magnitude.  When  u>2  increases 
through  the  value  given  by  (a),  a  pitchfork  bifurcation  occurs  at  (^1,^2)  =  (0,0)  wherein  this 
equilibrium  is  transformed  from  being  an  index  zero  critical  point  of  L  to  an  index  one  critical 
point.  In  the  process,  two  index  zero  critical  points  are  spawned.  When  u>2  increases  through 
the  value  given  by  (b)  (resp.  (c)).  the  equilibrium  (vr, 0)  (resp.  (0,tt))  is  transformed  from 
having  index  one  to  index  two.  Two  new  eqilibria  of  index  one  are  spawned  at  this  bifurcation. 
Finally  when  u>2  increases  through  the  value  given  in  (d).  the  equilibirum  (0,0)  goes  from 
index  one  to  index  two.  giving  rise  to  two  new  index  one  equilibria  in  the  process. 


(ii)  The  generic  bifurcations  in  a  gradient  system  such  as  this  depending  on  three  param¬ 
eters  r,  a,  a>2  involve  saddle  nodes.  (See  e.g.  Guckenheimer  and  Holmes.  1983.)  Each  time 
a  bifurcation  occurs  for  particular  parameter  values  and  (^1,^2)  =  (^i  > )»  simple  physical 
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considerations  dictate  that  a  symmetric  bifurcation  occurs  at  the  same  parameter  values  and 

(V’i.V’z)  =  (— -v>2)- 


As  u )2  —*  oo.  it  follows  as  in  lemma  2.2  that  solutions  to  (3.1)  approach  solutions  to 

(sin  V1!  +  txr  sin  rj> 2 )  cos  =  0 


(sin  V»i  +  r  sin  )  cos  ^2  =  0.  (3.2) 

The  solutions  to  this  system,  together  with  the  index  of  the  Hessian  evaluated  at  each  solution 
are  given  in  Table  3.1.  The  conclusion  of  the  theorem  is  a  straightforward  consequence  of 
these  facts. 


Table  3.1(a):  The  16  solutions  to  (S.t)  in  case  ar  <  1. 
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1  Cl 

Si 

C2 

*2 

Index 

0 

1 

0 

1 

0 

0 

-1 

0 

1 

0 

0 

1 

0 

-1 

1 

0 

-1 

0 

-1 

1 

0 

1 

v'i  -  iA2 

-l/r 

1 

0 

1 

i 

i 

1 

M 

I 

to 

-l/r 

1 

0 

-1 

,/l  -  l/r2 

l/r 

1 

0 

-1  v/1  -  l/r2 

1 A 

1 

1 

0 

1 

0 

2 

1 

0 

-1 

0 

2 

-1 

0 

1 

0 

2 

-1 

0 

-1 

0 

2 

Table  3.1(b):  The  11  solutions  to  (S.t)  in  case  ar  >  1. 


We  conclude  by  considering  the  special  case  r  =  3/2.  Using  MACSYMA  we  are  able  to 
determine  curves  in  the  (a,w)-parameter  space  where  bifurcations  occur.  These  are  depicted 
in  Figure  3.1  where  we  have  normalized  l  to  be  1.  Table  3.2  shows  how  the  distribution  of 
critical  points  of  each  index  is  related  to  the  total  number  of  solutions  to  (3.1). 


Figure  3.1:  Bifurcation  locus  of  equations  (S.l)  in  the  case  r  =  1.5.  The  top  curve  is 
cusped  as  approximately  a  =  .33. 
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t 

Ci 

Ci 

Ci  Ci 

Ci 

0 

1 

2 

2  2 

2 

1 

2 

3 

4  5 

6 

t 

1 

1 

2  3 

4 

total  4 

6 

8  10 

12 

Table  3.2(a): 

When  ar  >  1  the  number  of  critical  points  of  L  (=  solutions  to  (3.1)) 

increases  with 

|a>|  from  4  to  12. 

The  table 

gives  the  number 

of  critical  points  of  each  index. 

i 

Ci  Ci 

Ci 

Ci 

Ci  Ci 

0 

1  2 

2 

2 

t  4 

1 

2  3 

4 

5 

6  8 

2 

1  1 

2 

3 

4  4 

total 

4  6 

8 

10 

12  16 

Table  3.2(b): 

When  ar  <  1  the  number  of  critical  points  of  L  also  depends  on  |u>|.  This 

is  no  longer  a 

monotonic  dependence, 

however,  as  indicated  by  the  cusped  bifurcation  lo- 

cus  illustrated 

in  Fig.  3.1.  The  table 

gives  the  number  of  critical  points  of  each  index 

cooresponding  to  each  possible  number  of  solutions  to  (3.1). 


This  example  and  the  approach  we  have  used  to  study  it  lead  us  to  several  observations 
about  the  dynamics  of  more  general  chains.  We  remark  that  for  the  two  link  chains  under 
investigation  there  appear  to  be  a  less  interesting  set  of  bifurcations  when  r  <  1.  (The 
possibility  of  16  equilibria  as  u>  — *  oo  is  ruled  out.)  For  r  =  1.  a  bifurcation  occurs  at  u>  =  oo, 
as  illustrated  in  Table  3.1  where  several  solutions  to  the  equilibrium  equations  coalesce.  In 
general,  at  u  =  oo  we  may  in  principle  explicitly  display  the  4n  solutions  to  (2.2).  In  proving 
Theorem  3.2  in  Baillieul  (1987).  it  was  established  that  for  *  =  oc  and  generic  choices  of  the 
mass  and  link  length  parameters,  there  are  (£)  -2”  solutions  to  (2.2)  with  x,  ±  0  for  exactly  k 
values  of  t  between  1  and  n.  Certain  of  these  solutions  must  coalesce  when  all  link  lengths,  ii. 
and  mass  parameters,  m;,  are  equal.  To  see  this,  observe  that  when  =  tj  =  •  •  •  =  ln  =  l 
and  mi  =  m2  =  •  •  •  =  mn  =  m.  the  system  (2.2)  at  u  =  oo  is  equivalent  to: 

(nyi  +  (n  —  l)j/2  +  •  •  •  +  =  0 

((n  —  l)yi  +  (n  -  l)y2  + - 1-  y«)*2  =  0 


(yi  +  y2  +  — h  yn)*n  —  o 
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*,?+y,?-l  =  0  (t  =  1, . . .  ,n). 

Consider  solutions  of  the  form  xi  =  0  and  x*  ^  0  for  *  =  2, ...,ra.  For  this  case,  yi  =  ±1 
and  the  above  system  of  equations  reduces  to 

(n  -  l)j/2  +{n-  2)ys  +  •  •  •  +  yn  =  T(«  -  1) 

(n  —  2)yj  +  (n  —  2)ys  + - 1-  yn  =  T(n  ~  2) 


V2  +  1/3  H - l-yn  =  Tl- 

Subtracting  the  second  equation  from  the  first,  we  obtain  ya  =  =pl  and  X2  =  0.  showing  that 
this  solution  has  coalesced  with  the  solution  in  which  we  a  priori  stipulate  z2  =  0. 

That  there  are  fewer  than  the  maximum  possible  number  of  distinct  equilibria  when  we 
have  identical  link  lengths  and  joint  masses  suggests  that  there  may  a'so  be  fewer  stable 
equilibria.  We  may  obtain  insight  into  the  situation  for  large  values  of  u  from  the  following. 

Theorem  3.2:  For  an  n-link  chain  of  the  form  under  discussion  with  equilibrium  dy¬ 
namics  given  by  (2.2),  the  only  stable  equilibria  at  w  =  oo  are  those  for  which  each  is 
equal  either  to  y  or 

The  proof  of  this  theorem  will  not  be  given  here,  but  the  key  idea  is  that  if  some  ipi's 
take  values  other  than  ±|.  the  linkage  will  be  equivalent  in  equilibrium  to  a  chain  with  fewer 
links  which  is  in  an  unstable  equilibrium. 


Theorem  3.3:  For  an  n-link  chain  (as  depicted  in  Fig.  2.1)  with  equal  link  lengths 
=  ...  =  £„=£)  and  equal  joint  masses  (mj  =  •••  =  m„  —  m),  the  only  stable  equilibria 
at  w  =  oo  are  those  for  which  all  rpi  =  j  (or  symmetrically  all  ipi  =  —  y . 

Proof:  By  theorem  3.2  we  know  that  at  w  —  oo  all  V’*  at®  if-  Thus  we  may  assume 
without  loss  of  generality  that  for  some  integer  k  <  n.  V’n-fc  =  f  and  V’n-k+i  =  i’n-k+2  = 
■  ■ .  —  ij>n  =  —  At  u>  =  oo.  when  each  V’*  =  if .  the  Hessian 


92L 

dip2 


—  diag[(nsi  +  (n  —  1)<2  +  *  •  •  +  sn)Tnt2  Si , 


((n  -  l)aj  +  (n  -  1)j2  H - 1-  »n)m£2s2, 

■  •  • ,  (aj  -f  *2  +  •  •  •  +  ai»)tn/aan]- 
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For  stable  equilibria,  each  of  these  diagonal  entries  will  be  positive.  Hence  we  have 


(k  4-  l)(si  +  •  •  •  +  *n-k)  +  ksn-k-i  + - b  an  >  0 


while 


&(*!  +  ■  •  ■  +  an-fc)  +  kan-k-l  +  •  ■  •  +  «n  <  o. 


Let  a  =  H - b  an_*.  These  inequalities  may  be  rewritten  as 


( k  +  l)a  > 


k(k  +  1) 


ka  < 

which,  taken  together,  are  equivalent  to 


2 

k(k  + 1) 


k  k  + 1 

-  <  a  <  -  — . 

2  2 

But  this  combined  inequality  can  never  be  satisfied  for  integral  values  of  a.  Because  of  this, 
we  conclude  that  the  only  stable  configurations  occur  when  either  all  V’.  =  f  or  all  V**  =  —  f 
for  (x  —  1, . . . , n).  This  proves  the  theorem. 


4.  CONCLUSION 

Theorem  3.3  suggests  that  for  all  sufficiently  large  values  of  slew  rates  chains  with 
equal  link  lengths  and  joint  masses  will  have  only  one  (symmetric  pair)  stable  equilibrium. 
This  will  occur  with  the  chain  being  flared  out  with  all  angles  fa  having  the  same  sign.  As 
was  indicated  in  the  special  case  n  =  2,  in  order  for  there  to  be  a  richer  set  of  dynamically 
stable  equilibria,  uneven  mass  sizes  and  unequal  link  lengths  are  called  for.  While  the  models 
we  have  studied  are  idealizations  of  physical  chains,  our  laboratory  experiments  with  the 
controlled  slewing  of  actual  chains  seem  to  validate  the  theoretical  predictions. 
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